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ABSTRACT 

We select a sample of 70378 E/SO (early-type) galaxies at < z < 0.36 from the Sloan 
Digital Sky Survey, excluding disk and star-forming galaxies. We estimate g and r 
magnitudes in the observer- and rest-frames directly from the SDSS DR6 spectra; 
this provides an object-by-object estimate of the fc-correction. Observer-frame colours 
can also be estimated from the imaging photometry. However, in this case, rest-frame 
colours, and hence fc-corrections, must be inferred from fitting to data in other band- 
passes as well as to stellar population synthesis models. There are small (a few 0.01 
magnitudes) discrepancies between the spectra and imaging photometry, particularly 
for galaxies with with low signal/noise ratios in the spectra and more negative eclass 
spectral classifications. We correct for these, and then use the fc-corrections from the 
spectra to study the evolution of the rest-frame colour-magnitude (CMR) and colour-cr 
(CcrR) relations. 

Both the CMR and CcrR relations evolve blueward with increasing redshift, ap- 
proximately in agreement with passive evolution models with age ~ 12 Gyr. The rate 
of evolution is sensitive to the fc-corrections. Using fc-corrections from the CWW tem- 
plate spectrum, rather than the observed SDSS spectra, gives a CMR with much more 
evolution; fc-corrections from Blanton & Roweis, which are based on fitting Bruzual & 
Chariot models with solar abundance ratios to the ugriz colours, give less/no evolu- 
tion. However, the slope and zero-point of CMR relation depends on whether colours 
were defined in fixed physical or angular apertures, a consequence of the fact that the 
centers of these objects tend to be redder: the relation is steeper for fixed angular aper- 
tures. These colour gradients must be accounted for when estimating CMR evolution 
with fixed angular apertures. On the other hand, although the zero-point of the CctR 
depends on the aperture in which the colour was defined, the slope does not, suggest- 
ing that colour gradients are correlated with residuals from the a — Mr relation. Since 
these residuals are age indicators, our findings suggest that colour gradients depend 
on the age of the stellar population. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters 
- galaxies: evolution 



1 INTRODUCTION 

It has long been known (e.g. Visvanathan and Sandage 1977) 
that the more luminous early-type (elliptical and SO) galax- 
ies tend to have redder colours. Expressed in magnitudes, 
the colour-magnitude relation (CMR) is approximately lin- 
ear. The form of the CMR is most readily observed in a rich 
galaxy cluster, which provides a very wide luminosity range 
of E/SOs at a single redshift. 

The CMR provides significant information about galaxy 
evolution. Firstly, from the tightness (scatter < 0.05 mag) 
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of the CMR for the E/SOs in the Coma and Virgo clusters. 
Bower, Lucey and Ellis (1992) concluded that these galaxies 
had formed at z > 2. Secondly, the slope of the CMR ap- 
peared to derive primarily from a correlation of metallicity 
with galaxy mass (e.g. Kodama and Aromoto 1997; Vazdekis 
et al. 2001). Higher metallicity results in redder colours for 
a given age. The observed trend would reflect the greater 
ability of more massive E/SOs to retain metal-enriched gas 
during their formative starbursts. An increase in stellar age 
with galaxy mass may play a secondary role. 

A 'red sequence' of E/SO galaxies is seen in rich clusters 
out to at least z = 1.27, with little or no change in the CMR 
slope (Blakeslee et al. 2003; Mei et al. 2006, 2009). The red 
sequence appears to have first formed at z ^ 2 (Labbe et al. 
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2007). For a passively aging stellar population there would 
be a blueward shift of the CMR with increasing lookback 
time (redshift); Holden et al. (2004) and Taylor et al. (2008) 
find some evidence of this colour evolution. 

In this paper we aim to obtain an accurate measurement 
of the evolution of the CMR to 2; = 0.36. However, to de- 
termine the CMR over a wide redshift range, but in a fixed 
rest-frame colour (e.g. g — r), we need to apply accurate 
fc-corrections to convert the observed photometric magni- 
tudes back to the galaxy rest-frames. Galaxy fc-corrcctions 
are commonly derived from modelled spectra, but in reality 
each galaxy has its own individual /c-correction, and the un- 
certainties in modelled fc-corrcctions (even a few 0.01 mag) 
may well be comparable to the small evolution that we aim 
to study in the CMR. Thus, it would be much preferable 
to derive either fe-corrections, or rest-frame colours, directly 
from the galaxy spectra, but for good statistics on the CMR 
a database of many thousand high quality, flux calibrated 
spectra are required. 

The Sloan Digital Sky Survey (SDSS, Stoughton et al. 
2002) had as one of its stated aims to revolutionize the study 
of galaxies at moderate redshifts, by obtaining hundreds of 
thousands of spectra together with five-band (ugriz) pho- 
tometry. Bernardi et al. (2003a) carried out an initial study 
of the early-type galaxies in the SDSS, selecting 9000 E/SOs 
(from the ~ 65000 galaxy spectra available at that time) 
using a combination of morphological and spectral criteria. 
Rest-frame colours and magnitudes wore estimated, using 
fc-corrections based on (i) an evolving Bruzual & Chariot 
(2003, hereafter BC03) model with a single burst of star- 
formation 9 Gyr ago, with solar abundances and a Kroupa 
(2001) IMF, (ii) an observed elliptical template spectrum 
from Coleman, Wu and Weedman (1980, hereafter CWW). 

For these E/SO galaxies Bernardi et al. (2003b) esti- 
mated evolution rates of I^{Mg) = — 1.15z, A(Mr) = —0.85a 
and A(Mi) = -0.752. Bernardi et al. (2003c) derived a 
CMR using a k-correction from a composite of the above 
two models, Kh^^wi^) = Klwwi^)+K'Sc " ^Sc"'"- The 
rest-frame colour [g — r)rf correlated with luminosity and 
velocity dispersion as (-0.025 ± 0.003)M^ and (r°-^®=^° °^. 
Metal abundances in the E/SOs correlated with velocity dis- 
persion e.g. as Mg2 cx a^-^^^o-^'^ and (Fe) oc o-"■"±°■°^ and 
the mean stellar age was estimated as 9 Gyr. 

Bernardi et al. (2005) re-examined the CMR, using 
39320 early-type galaxies from the SDSS, and showed that 
it was the result of more fundamental correlations between 
colour and velocity dispersion a, and between a and lumi- 
nosity. 

Gallazzi et al (2006), using a similar sample of SDSS 
E/SOs, found a steep correlation of both mctallicity and 
Mg/Fe ratio with stellar mass, as traced by either lumi- 
nosity or velocity dispersion (e.g. Mg2 oc (j^-'^^^^-'^'^y They 
explained the CMR as a sequence in stellar mass, with some 
increase in mean age. Their CMR exhibited only a min- 
imal redward shift (0.006 ± 0.003 mag in p — r) between 
less and more dense clustering environments. Jimenez ct al. 
(2007) concluded from the small scatter in the CMR that 
the E/SOs formed most of their stars aX z > 2, and from 
the super-solar Mg/Fe abundance ratios ('a-cnhancemcnt') 
of the more massive E/SOs that their star-formation took 
place in a short burst (~ 1 Gyr). 

Continued and ongoing SDSS observations have greatly 



enlarged the database of spectra, enabling us to analyse the 
CMR of a larger sample of E/SOs than previously. In Section 
2 we describe the data and sample selection, and in Section 
3, the derivation of the fe-corrections from the spectra. In 
Section 4 we present our CMR and compare it with CMRs 
obtained using other apertures and fe-correction prescrip- 
tions. We present a similar analysis of the CcrR in Section 
5. We summarize and discuss our results in Section 6. 

SDSS magnitudes are given in the AB system where 
rriAB = -48.60 - 2.5 log F,. and ruAB = is 3631 .ly. We 
assume throughout a spatially flat cosmology with Ho = 70 
km s~^Mpc~^, Qm = 0.27 and Qa = 0.73, giving the age of 
the Universe as 13.88 Gyr. 



2 SDSS DATA 

2.1 Selecting the E/SO sample 

Our sample is selected from Data Release 4 of the Sloan 
Digital Sky Survey (Adelman- McCarthy et al. 2006), which 
covered 4783 deg^ of sky with spectra for a total of 
673280 sources. However, our analysis uses more recent ver- 
sions of the spectra from SDSS Data Release 6 (Adelman- 
McCarthy et al. 2008), calibrated at Princeton University 
(spectro.princeton.edu). These reductions include measure- 
ments of the line-of-signt velocity dispersion (cr) through the 

3 arcsec diameter aperture. 

This dataset contained 367471 galaxies of all types (but 
excluding stars, QSOs and unclassifiable sources) in the 
range 14.5 < r < 17.77 (Petrosian AB magnitudes) with 
flux/wavelength calibrated optical spectra (approx. 3800- 
92OOA with resolution 1800-2100) and ugriz photometry. 
The SDSS image processing software provides several sets 
of ugriz magnitudes for each object. The photometric sys- 
tems of interest here arc: 

(i) Model magnitudes: the total magnitude corresponding 
to the best-flt profile, which is for these galaxies a de Vau- 
couleurs (1948) profile, with an effective radius (rdev) fitted 
in the r band. The radius rdeV is kept fixed for the mag- 
nitude measurement in the other passbands, so as to give 
accurate colours even for galaxies with a strong colour gra- 
dient. 

(ii) Fiber magnitudes: measured in fixed 3 arcsec apertures 
on the imaging data after smoothing to a uniform 2 arc- 
sec resolution - intended to correspond to the spectroscopic 
fiber aperture. 

(iii) Petrosian magnitudes. The Petrosian (1976) radius rpet 
is the angular radius where the mean surface brightness av- 
eraged over the aperture r < rpet is 5 times higher than 
the mean surface brightness on the annulus r — rpet- The 
Petrosian flux/magnitude is the total flux observed within 
twice this radius. 

We apply strict selection criteria in order to extract a clean 
sample of E/SOs with good spectra, and exclude spirals, 
galaxies with disk components, and star-forming galaxies. 
The most important selection criteria are that: 
(i) fracdcv(r): the fraction of total flux contributed by the 
de Vaucouleurs component when the galaxy profile is fitted 
by the sum of de Vaucouleurs and exponential profiles. We 
include only non-disk galaxies with fracdev of unity for both 
the g and the r band. 
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(ii) Eclass, a ID classification of spectral type, obtained 
from a principal component analysis (Yip et al. 2004) - lower 
(negative) values corresponding to absorption-line galaxies 
with old stellar populations, positive values indicating star- 
forming galaxies. We include only galaxies with eclass < 0. 
In addition we require that: 

(iii) Concentration, defined as the ratio of the Petrosian ra- 
dius rpet, 90 (the angular radius containing 90% of the Pet- 
rosian flux) to the smaller radius rpet,50 (containing 50%), 
measured in the i-band, must be > 2.4 to reduce possible 
contamination of later-type galaxies (this selection excludes 
only a few objects). 

(iv) The r band model magnitude, corrected for Galactic 
reddening, must be 14.5 < r < 17.5 (this excludes some of 
the faintest objects). 

(v) The velocity dispersion (a) as measured from the SDSS 
spectrum must be at least 60 km s^^, but no more than 450 
km s~^ (the upper limit to exclude superimposed galaxies, 
as no single galaxy in the SDSS appears to have a > 450 
km s~^; e.g Bernardi et al. 2008; Salviander et al. 2008) 

(vi) Redshift must be 0.005 < z < 0.36. The reason for the 
upper limit is that, at z > 0.36, the rest-frame r-band be- 
gins to shift significantly outside of the 3800-9200A range of 
the SDSS spectra, so we cannot obtain an accurate r-band 
k-correction. 

The selected galaxies, satisfying all the above, number 
70378. 



2.2 Corrections applied to SDSS parameters 

We apply a number of, generally very small, corrections to 
the a, radii and model magnitudes given by the SDSS cat- 
alogs. 

The a measured for early-type galaxies decreases (very 
slowly) with distance from the centre, and hence the SDSS 
a will have a dependency on the ratio of the spectro- 
graph aperture to the galaxy's effective radius. Jorgensen 
et al. (1995) estimated that a cx {rap/rdev)~'''^'^ , where 
TdeV is the de Vaucouleurs model radius in the r-band. As 
in Bernardi et al. (2003a), we correct all the a measure- 
ments from the SDSS spectra to apertures of Vdev/S. As 
rap ~ 1.5 arcsec for the SDSS, and VdeV is in arcsec, this 
gives acorr = t^SDSsirdeV /l'2-0)~°'°* ■ 

We correct the de Vaucouleurs model-fit radii, which 
are semi-major axes, to effective radii (reff) by multiplying 
by the square root of the axis ratio, r^// = rdev\/b/a. 

Furthermore, there is evidence of systematic errors in 
the sky subtraction in the SDSS reductions (e.g. Bernardi et 
al. 2007; Lauer et al. 2007; Hyde & Bernardi 2009). These 
caused the size and total flux of large, extended objects to 
be underestimated. By comparing SDSS catalogs and more 
accurate reductions, Hyde & Bernardi (2009) estimated cor- 
rections for the SDSS effective radii. These were positive, a 
steep function of r^ff, and zero for Veff < 2.0 arcsec. Cor- 
rections were also estimated for the de Vaucouleurs model 
fit magnitudes, these being also a function of re//, but zero 
for reff < 1.5 arcsec. 

We correct all our re// radii as in Hyde & Bernardi 
equation 4. The model magnitudes, for these galaxies with 
fracdcv ~ 1-0, are the same as de Vaucouleurs model-fit 
magnitudes, except that the radius of the r-band fit is used 



for all passbands. Hence we correct our model magnitudes in 
all (ugriz) passbands as in equation 3 of Hyde & Bernardi 
(2009) but always using the r-band r^ff. This correction will 
therefore cause no change in the model-magnitude colours. 
For the fiber magnitudes, no sky-subtraction correction is 
applied, as these small aperture magnitudes correspond ap- 
proximately to a model fit with re// — 1.5 arcsec, where the 
this correction falls to zero. 



3 DERIVING fs:-CORRECTIONS FROM 
SPECTRA 

3.1 Method 

We calculate g and r-band fc-corrections directly from 
the fiux-calibrated spectra, corrected for atmospheric and 
Calactic reddening. To determine the k-correction for a 
galaxy at redshift z, we first integrate the galaxy spectrum 
F{\) over the filter response function. In the g-band, g{X), 
the apparent magnitude Qspec, in the AB system, is (Hogg 
2002; Bruzual & Chariot 2003) 



gspec = -2.5 loglO 



/ F{\)g{X)\^ dlnX 
J C{X)g{\)X^dlnX 



(1) 



where C(A) = 3631 (c/A^) x 10"^^ erg s"^ cm"^ is the spec- 
trum of a source with = 3631.0 Jy at all frequencies, 
defining rriAB ~ 0. In practice, this integration is performed 
by summing over the SDSS pixels which are given as Fx at 
intervals of A (log A) — 0.0001. If we sum over the pixels, 
then 



gspec = -2.5 loglO 



J^iF.Wi X\lc 



3631 X 10-23 erg s-^cm" 



(2) 



where the summation is over all pixels (3800-9200 A), and 
we have defined Wi = gi/Tngi. 

We then calculate the magnitude gospec which would 
be measured in the p-band if the spectrum is de-redshifted 
by a factor 1 + z, back to 2 = 0. To conserve flux, this 
'compression' increases the F\ normalization by a factor 1-1- 
z, making 

gospee -2.5 loglO 3631 ^ 10-23 erg S-W2 ' 

We define Vspec and rospec similarly, so the ^-corrections are 



kg — QOspec Qspec and kr — ^Ospi 



ec ' spec- 



(4) 



3.2 Spectral and photometric calibrations 

Magnitudes derived from the spectra are compared with the 
(reddening-subtracted) fiber (i.e. 3 arcsec aperture) mag- 
nitudes from the DR6 SDSS catalog. It is expected that 
these would be offset by ruapec — rrifib — —0.31, because 
the spectra are flux-calibrated using the PSF magnitudes of 
the spectrophotometric standard stars (Adelman-McCarthy 
et al. 2008), rather than the fiber magnitudes (which are 
about 0.31 mags fainter). 

Figure [T] shows a histogram of nispec — mfu,, in the 
g, r and i bands, for our E/SO sample, in a few redshift 
bins. At 2 < 0.15, we find {vspec — rfit) = —0.299 and 
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Figure 1. Histogram of the difference between the spectra- 
derived and fiber (3 arcsec aperture) magnitude of SDSS early- 
type galaxies, plotted for the g, r and i bands and divided into 3 
redshift intervals. 



{ispec — ifib) = —0.308, but (gspec — gfib) — —0.259. This dis- 
crepancy means that the g — r extracted from the spectra are 
on average 0.04 mag redder than the fiber magnitude g — r. 
The calibration of the g — r magnitudes from SDSS data is 
claimed to be accurate to ~ 0.02 mag (Adelman-McCarthy 
et al. 2008), but it is conceivable the errors from the imaging 
and spectroscopic data might be in opposite directions from 
the true calibration, producing a 0.03-0.04 mag difference 
in their g — r. 

In addition, Figure [2] shows that the difference between 
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Figure 2. The mean difference between the spectra-derived and 
fiber (3 arcsec aperture) magnitudes of SDSS early-type galaxies, 
as a function of redshift with the galaxies divided into four in- 
tervals of the seeing during the spectroscopic observation. Shown 
for g (top) and r (bottom) bands. 



the flux recorded in the spectra and that on the imaging 
data is affected slightly by the atmospheric seeing during the 
spectroscopic observation. Imaging data are taken in good 
seeing, but 51% of the spectra were taken in r-band seeing 
of > 1.8 arcsec. We estimate each 1 arcsec increase in seeing 
produces a flux deflcit A{mspec — rnjib) of 0.021 mag in g 
and 0.012 mag in r. Typical seeing could then redden the 
spectra by A((; — r) ~ 0.01 mag, but does not account for 
all the difference in spectra and fiber g — r. The relatively 
small effect of seeing on the spectrophotometry should not 
significantly affect our estimates of the CMR evolution as 
it does not appear to increase with redshift, or to fainter 
fluxes. 

However, Figures [1] and [2] also show that the mean 
gspec — gfib increases with redshift, as does the scatter in 
this quantity. This is not related to seeing but might be 
be a separate problem with the SDSS spectrophotometry. 
Further investigation of the gspec — gfib offset found it was 
correlated with the eclass spectral type classiflcation pa- 
rameter, being much greater for redder galaxies with more 
negative eclass (Figure [3]). Note the E/SO sample all have 
eclass < with a mean of -0.1275. In the red-band the mean 
r spec — r fib is also correlated with eclass, although the effect 
is smaller and does not increase strongly with redshift. The 
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Figure 4. Spectra of four representative z ~ 0.3 E/SOs with sim- 
ilar magnitudes and a range of eclass, showing the difference 
at the blue end, with the g-band response function in the ob- 
server frame and in the galaxy rest-frame. Spectra are smoothed 
(slightly) with a = 2 pixels. 



dependence on eclass is obviously non-linear and might be 
better described as eclass^. 

To understand why (/spec might be affected by eclass, 
especially at higher redshifts, we examine the spectra of 
four representative E/SO galaxies of different eclass, all at 
z > 0.25. Figure |4] illustrates that a more negative eclass is 
associated with a lower flux in the rest-frame UV, which at 
these redshifts is observed in the p-band. The spectra are 
particularly noisy at the blue end of the observed p-band, 
near 4000 A. 

Flux variations from noise at the blue end of the 
spectrum would produce correlated variations in Qspec and 
eclass. But the trend in {gspec — gfib) with redshift implies 
that this scatter is not symmetric, but skewed. As a result, 
for spectra with very low signal/noise at these wavelengths, 
the integrated pspec flux is, on average, underestimated by 
a few times 0.01 mags. 



3.3 Empirically improved calibration 

Although the reason for this discrepancy in the SDSS spec- 
trophotometry is not known, we can try to correct for it as 
follows. We begin by assuming that the offset between the 
spectra-derived and fiber magnitude is not directly sensitive 
to redshift or wavelength, but depends only on: 



(i) the signal/noise in the spectrum, averaged across the 
passband, 

(ii) the tilt of the spectrum across the band, represented by 
eclass , 

(iii) a cross-term (eclass^ x signal/noise). 

For each galaxy we calculate directly from the SDSS spec- 
trum and the associated noise spectrum, the root-mean- 
square signal/noise per pixel in 4 wavelength intervals, the 
observer-frame and the rest- frame g and r bands. Figure [S] 
shows {gspec — gfib) as a function of rms signal/noise, for 
four bins in eclass. The mean trend is well described by 

^^o-,n, 0.09186 0.505216 

{gspec - gfib) ~ -0.28731 \ — 



-0.32549y - 7.2538 



y 



(5) 



where x is the rms signal/noise and y is eclass^. The rms 
residual is 0.060 magnitudes. The smooth curve in the Figure 
shows this expression when eclass= —0.1275. 

We define our signal/noise correction function to be the 
final four terms in equation ((5)1 , 

Before proceeding further, we illustrate the effective- 
ness of this correction in removing the differences between 
the spectra and image photometry, in particular the redshift 
trend in Figure|3] Figure[6]is similar to Figure|3l except that 
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Figure 5. The mean gspec — 9 fib of early-type galaxies from the 
SDSS photometry, as a function of root-mean-square signal/noise 
pixel~^in the g-band of the spectrum. Galaxies are divided into 
4 intervals of eclass. A best-fit function of signal/noise (x) for 
the sample mean eclass of -0.1275 is also plotted. 



now the signal/noise corrections (based on equation[S| have 
been subtracted from both gspec and Vspec- Almost all the 
redshift and eclass dependence has now been removed from 
gspec — g f ib, although there may be slight residuals especially 
at 2 > 0.25. Similarly, by using the rms signal/noise in the 
rest-frame g and r, we obtain the corrections for gspeco and 
rspeco, and then a set of 'corrected' k-corrections for each 
galaxy from the gspec — gspecO and rspec — VspecO- 



3.4 Comparison with models 

We now compare our estimated fc-corrections with those de- 
rived from models which represent the range of expected 
ages and star-formation histories for E/SOs (e.g. Jimenez 
et al. 2007). We do this by using the Bruzual & Chariot 
(2003) spectrophotometric models, with the (recommended) 
Padova 1994 stellar evolution tracks, the Chabrier (2003) 
IMF and fixed solar metallicity. 

In our Model 1, star-formation (SF) begins 12.5 Gyr 
ago {z — 4.53), continues at a constant rate for 1 Gyr 
until 11.5 Gyr ago {z = 2.83) and then ceases. In Mod- 
els 2, 3 and 4, SF begins 12.0 Gyr ago at {z = 3.49) and 
the SF rate decreases exponentially with timescales t — 1, 

1.5 and 2 Gyr. The models include no internal dust ex- 
tinction after the cessation of SF. Over < z < 0.3, 
models (1,2,3,4) predict mean rates of luminosity evolution 
A{Mr) = -(1.11,1.20,1.37,1.52)2 and rest-frame colour 
evolution A{g - r) = -(0.124, 0.171, 0.265, 0.366)2. 

We can also estimate the evolution of the luminosity 
function characteristic magnitude M* directly from the data 
using a V/Vmax method. We sum V/Vmax with the rest- 
frame absolute magnitudes adjusted faintward by a linear 
function of redshift +Qz, which is a reasonable approx- 
imation at these moderate redshifts. The mean V/Vmax 
matches the no-evolution value of 0.50 for Q = 0.86, thus 
giving an estimate of the luminosity function evolution as 
A{Mr) — — 0.862. This agrees with the —0.852 estimate of 
Bernardi et al. (2003b) for SDSS early-types and the -0.902 
of Hyde & Bernardi (2009) . The difference from the slightly 
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Figure 6. The mean magnitude difference (a) gspec — g fiber (b) 
rspec — T^odel of SDSS early-type galaxies, as a function of red- 
shift, divided into 5 intervals of eclass, and using the 'corrected' 
gspec and rspec magnitudes. 



stronger evolution in the BC03 models could be explained 
by a small amount of merging. 

Figure [7] shows the mean observer-frame g — r colours 
(spectra-derived, model and fiber) of our E/SO sample as 
a function of redshift, compared to the BC03 models. No 
single model fits over the whole redshift range, because 
with increasing redshift, more intrinsically luminous galax- 
ies dominate and these (as we shall see) are significantly 
redder. The uncorrected spectra-derived colours gradually 
cross from those of BC03 model 4 at 2 ~ to the redder 
model 1 at 2 > 0.3, while the mean fiber and model mag- 
nitude colours are at all 2 a few 0.01 mag bluer, and closer 
to model 2 at 2 > 0.3. The signal/noise corrected colours 
are intermediate, being similar to the uncorrected colours 
at low redshift but meeting model 2 at 2 = 0.3. The mean 
model colours are slightly (< 0.02 mag) bluer than the fiber 
colours, especially at 2 < 0.1, which may be the result of ra- 
dial colour gradients. 

Figure [8] shows the spectra-derived k-corrections of the 
E/SOs - from the uncorrected and corrected magnitudes 
- averaged in A{z) = 0.025 intervals, and compared with 
the four BC03 models and CWW template spectra. The 
k-corrections are tabulated in Tables 1&2. Note that at red- 
shift 2, the BC03 model k-correction depicted is that cal- 
culated on the model spectrum at the lookback time t{z), 
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Figure 7. Mean observer-frame g — r colours of the E/SO galaxies 
as a function of redshift, calculated from the spectra — before 
and after the signal/noise 'correction' - and from the fiber and 
the model magnitudes. These g — r colours are compared with 
four evolving BC03 models for early- type galaxies (see text). For 
clarity a linear trend 2.52; is subtracted in all cases. 



and not the k-correction of the z = spectrum, i.e. we plot 
the quantity Mrf{z) — Mev{z), and not 'fc' from the 'mag- 
nitude_FN' file (see BC03). 

The signal/noise correction reduces the (7-band k- 
correction (kg) at higher redshifts, by ~ 0.05 magnitudes 
at z = 0.3, but the reduction in kr remains very small at 
< 0.01 magnitudes. The mean fc-corrections from the spectra 
are generally in the range of the BCQ3 model predictions, 
nearest to Model 3 or 4 at lower redshift and the redder 
Models 1 (uncorrected) or 2 (corrected) at 2 > 0.2. 



4 COLOUR-MAGNITUDE RELATION (CMR) 

We now study the colour-magnitude relation and its evolu- 
tion. In all cases, the absolute magnitude Mr we use is 



Mr 



(6) 



where r is the model r magnitude from the photometry, Ar 
is the Galactic reddening in the r-band at the sky position 
of the galaxy (as given in the DR6 catalog), kr is the r- 
band k-correction, and dmod{z) the distance modulus at the 
spectroscopic redshift z. We discuss the effects of different 
choices for the fc-correction below. The final term Qz, with 
Q — 0.86 here, 'de-evolves' the absolute magnitudes to 2 = 
0, so as to isolate the effect of colour evolution from that of 
luminosity evolution. 

The rest-frame colour. 



{9 - r)rf = g - Ag - kg - {r - Ar ~ kr) 

= {g-r)-{Ag-Ar)-{kg-kr), 



(7) 



can be calculated in a number of different ways. If we use the 
spectra to estimate {go —ro) = (g — r) — {kg — kr), then these 
are colours within a fixed angular aperture. We compare this 
with the result of using g — r from the fiber photometry, 
but kg — kr from the spectra. We then show what happens 
if we determine the A;- correct ions directly from the fiber 
colours. 
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kg and kr, 

calculated from both uncorrected and corrected spectra magni- 
tudes, compared with four evolving BC03 models and also and 
the fc-correction from the CWW non-evolving template spectrum 
for E/SOs. 



To explore the dependence on aperture, we also use 
model g — r from the photometry, which gives colours within 
a physical scale which is related to the r— band half-light 
radius of each galaxy. In this case, if we determine the 
fc-correction from the spectra, then we have colours from 
one scale in the galaxy, but fc-corrections potentially from 
another. Therefore, we also show the result of deriving k- 
corrections from the model photometry, as was done for the 
fiber photometry. 

In all cases, we did not find evidence for evolution in 
CMR slope with redshift. On the basis of this, and previous 
studies finding no CMR slope evolution even to 2 ~ 1.3 (e.g. 
Mei et al. 2009), we fit (using IRAF surf it) 



(5 



Mr 



• ao + ai Mr + a2 2, 



(8) 



with separate linear dependencies on Mr and z and no cross- 
term. A few galaxies with outlying colours {g — r)rf < 0.5 or 
> 1.1 were excluded from the slope and evolution fits (0.55% 
of the total, leaving 69987 objects). 

4.1 CMR for a fixed angular scale 

Figure[9]shows the CMR from {g — r)rf colours directly mea- 
sured from the spectra, with and without the signal/noise 
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Figure 9. Colour-magnitude relation {g — r against Mr in the 
rest-frame) for SDSS E/SO galaxies in seven A(z) = 0.05 intervals 
of redshift (0.30 < z < 0.36 for the seventh), calculated using rest- 
frame g and r derived directly from the SDSS spectra, before (a) 
and after (b) applying the signal/noise correction (equation[5]and 
associated discussion). For clarity, points corresponding to fewer 
than 10 galaxies are not plotted. 



correction. The galaxies are divided into 7 redshift intervals 
to illustrate the evolution. The signal/noise correction has 
much less effect on the rest-frame gspec than on the observed- 
frame Qspec, and consequently it produces much less shift in 
the spectra-derived CMR than the g-hand k-correction. The 
colours of higher redshift galaxies are shifted bluewards, but 
by no more than 0.01 mag. The CMR, before and after cor- 
rection, is well described by equation (|8]) with 

ao = (0.2825 ± 0.0086) or (0.2426 ± 0.0083) 

ai = -(0.02486 ±0.00042) or - (0.02659 ± 0.00041) 

az = -(0.2368 ±0.0069) or - (0.2709 ± 0.0067) 

with an rms scatter of 0.0541 magnitudes (Figure |9}. Notice 
that the correction slightly increases the CMR evolution; it 
has little effect on the slope. 

On the other hand, the CMR appears slightly curved. 
This might cause a linear fit to mis-estimate the colour evo- 
lution. Fitting a quadratic 
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Figure 10. Same as previous figure, but now g — r is from ob- 
served fiber magnitudes, shifted to rest-frame using k-corrections 
derived from the spectra before (a) and after (b) applying the sig- 
nal/noise correction. 



yields 

ao = 

ai = 

a2 = 

as = 



(0.5119 ± 0.1066) or (0.1722 ± 0.1035) 
-(0.0985 ± 0.0099) or - (0.0650 ± 0.0096) 
-(0.00171 ± 0.00023) or - (0.00089 ± 0.0022) 
-(0.2227 ±0.0072) or - (0.2635 ± 0.0069) 



)rf 



■ ao + ai Mr + a2 Mr + as z, 



(9) 



This reduces the evolution only slightly. 

Figure [TOl shows the CMR from fiber (3 arcsec aper- 
ture) magnitudes, corrected to rest-frame using the k- 
corrections from the spectra. In this case, applying the sig- 
nal/noise based correction makes a big difference. We find 

ao = (0.2465 ± 0.0092) or (0.2685 ± 0.0094) 

ai = -(0.02517 ±0.00045) or - (0.02363 ± 0.00046) 

aa = -(0.3480 ±0.0074) or - (0.1999 ± 0.0075) 

before and after the correction. The uncorrected kg and kr 
yields very strong evolution; the correction reduces this sub- 
stantially (about 0.15z;). 



4.2 CMR for a 'fixed' physical scale 

Figure fTTI shows the CMR obtained using the observed model 
colours corrected to rest-frame with the spectra-derived k- 
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Figure 11. Rest-frame colour-magnitude relation in which g — r 
is based on model colours from the photometry and (signal/noise 
corrected) fc-corrections from the spectra. 
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Figure 12. Rest-frame colour-magnitude relation in which g — 
is based on de Vaucouieurs profile fits to the galaxies in g and 
and (signal/noise corrected) fc-corrections from the spectra. 



corrections. This relation is well fit by 



ao = (0.3543 ± 0.0091) or (0.3794 ± 0.0093) 

ai = -(0.01904 ±0.00045) or - (0.01735 ± 0.00046) 

aa = -(0.2242 ±0.0074) or - (0.0743 ± 0.0075) 

with an rms scatter of 0.0578 magnitudes. Notice that the 
rest-frame colours are ~ 0.04 mag bluer than before, and 
the evolution is much smaller when the signal/noise fc- 
corrections are used. In addition, the CMR is shallower and 
less curved. 

There is nothing physically significant about the size 
of the SDSS fiber, but the half-light radius may be a phys- 
ically meaningful scale - perhaps this is why it shows less 
curvature than the CMRs based on a fixed angular aperture. 
The fact that this CMR is bluer than that associated with 
a 3 arcsec aperture is due to colour gradients: the typical 
half-light radius of the objects in our sample is larger than 
3 arcsec, and E/SOs tend to be redder in the center. Many 
previous CMR studies e.g. van Dokkum et al (2000) and 
Mei et al. (2006, 2009) used colours measured in half-light 
apertures, so as to minimize colour gradient effects. The fact 
that the slope becomes shallower for the model colours may 
have implications for how gradients scale with luminosity, 
which we return to shortly. However, the fact that the evo- 
lution in this case is smaller than for a fixed angular aper- 
ture indicates that gradients have contributed to our other 
two estimates of the evolution of the CMR. (A fixed angu- 
lar scale collects light from a larger physical radius, where 
the colours are expected to be bluer even if there were no 
evolution, for the high redshift objects). 

Thus, there are many senses in which this CMR is 
the one with the most direct physical information. Unfor- 
tunately, it is based on fc-corrections which are computed 
for a different physical scale than that on which the colour 
itself is measured. The next subsection studies the magni- 
tude of this effect. 

Yet another measure of 'total' magnitudes and colours 
may be obtained by fitting de Vaucouieurs profiles indepen- 
dently in g and r, without constraining the effective radii to 
be the same. Figure [12] shows the CMR from this type of 
photometry (with the signal/noise correction). This can be 



cww - 

BC03:1 - 
BC03:2 
BC03:3 
BC03:4 - 
spectrum-derived ^ 
s-d [corrected] 
Blanton (fibre) - 
Blanton (model) - 




0.15 0.2 
redshift 

Figure 13. The difference kg — kr, with a linear —2.bz trend sub- 
tracted, for E/SO galaxies in A(z) = 0.025 bins, from the spectra 
(with and without the signal/noise correction), from Blanton & 
Roweis (2007) for model and fiber magnitudes), from the 4 BC03 
models, and from CWW. 



fitted with 

ao = (0.3455 ±0.0121) 
ai = -(0.01726 ±0.00060) 
aa = -(0.1558 ± 0.0098) 

and has a colour zero-point 0.03-0.04 mag bluer than the 
corresponding model magnitudes CMR. This again is due 
to colour gradients, which cause g-band effective radii to be 
systematically larger than those fitted in the r-band. How- 
ever, the CMR slope is not significantly different, which 
suggests the mean ratio of (observer-frame) g and r radii 
does not strongly depend on luminosity. The greater blue- 
ward evolution of the de Vaucouieurs CMR may mean that 
this ratio increases with redshift (from gradient evolution, 
or from steeper gradients at shorter rest-frame wavelengths, 
e.g. u — b). 
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Figure 14. Same as previous figure, but now sliowing the cor- 
rected spectra-derived k-correction for E/SO galaxies divided into 
4 intervals of absolute magnitude M^. 



4.3 CMR from other fc-corrections 

In this subsection, we consider the k-correction computed 
following Blanton & Roweis (2007), a program which fits the 
broad-band (ugriz) magnitudes of each individual galaxy 
with one of a set of BC03 model template spectra, from 
which the g and r-band k-corrections are then measured. 
We also compare with two simpler models in which the 
same k-corrections are assumed for all the E/SOs: firstly, 
our evolving BC03 model 2 (age 12 Gyr r = 1 Gyr) and 
secondly, the non-evolving for ellipticals from Coleman, Wu 
and Weedman (1980, hereafter CWW). Figure [13] shows the 
different kg — kr against redshift (Table 1&2 list the difi'erent 
fc-corrections used in the Figure). 

Firstly, note that k-corrections computed from model 
and fiber magnitudes differ (primarily at 2 < 0.15) by less 
than 0.01 mag. At 2; < 0.15 the spectra and Blanton- Roweis 
kg — kr are all below our BC Model 4, but the uncorrected 
kg — kr crosses all the models to reach Model 1 at 2 > 0.3. 

kr is lower at high redshift, and the 
- kr are lower still, and give the least 



The corrected kg 
Blanton-Roweis k„ 



that from the spectra, and so gives an overestimate of A{g~ 
r) and too-blue rest-frame colours. 

Neither the corrected or uncorrected kg — kr follow any 
of the BC03 models over the whole redshift range. This is 
to some extent because the mean luminosity increases with 
redshift. To examine the effect of this. Figure [TJ] shows mean 
kg — kg with the galaxies divided into four intervals of Mr- 
At a given redshift the more luminous galaxies do tend to 
have a slightly more positive kg — kr- However, the galaxy 
kg — kr within each of the four luminosity intervals still do 
not closely follow any of the BC03 models, tending to be 
lower at 0.05 < 2 < 0.2, indicating systematic differences 
at some wavelengths between the spectra and this set of 
models. 

We compared the spectroscopically-based CMRs with 
the CMRs from photometric colours fc-corrected using only 
model-based or template spectra. The top panel in Fig- 
ure [15] shows the effect of fitting Blanton & Roweis (2007) 
fc-corrections to fiber magnitudes, and using the observed 
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Figure 15. Colour-magnitude relation using fc-corrections from 
Blanton & Roweis (2007) fitted to ugriz fiber (top) and model 
magnitudes (bottom). 



fiber g-r. The CMR is best-fit with 

ao = 0.2852 ± 0.0080 

ai = -(0.0227 ±0.00039) 

aa = -(0.10541 ±0.0064); 

The slope is very similar to that in the bottom panel of 
Figure 1101 but with these rather than the spectra-based fc- 
corrections, the evolution is 0.092 weaker. A similar analysis 
of the model colour (Blanton & Roweis fc-corrections from 
the model ugriz magnitudes) is shown in the bottom panel of 
Figure [TSl Again the model-magnitude colours give a CMR 
with a shallower slope, 

ao = 0.3781 ± 0.0075 

ai = -(0.01761 ±0.00037) 

aa = -(0.0029 ±0.0060), 

but with these k-corrections the evolution is reduced to zero 
or even inverted! This may be an artifact of the fact that 
the BC03 models used by Blanton & Roweis assume solar 
abundance ratios, and this is unrealistic. 

Figure [16] shows the CMR from model colours and fc- 
corrections from BC03 Model 2 k-corrections. It is well fit 
by equation ([8]l with 

ao = 0.3241 ± 0.0094 

ai = -(0.01846 ±0.00046) 
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Figure 16. Colour-magnitude relation for E/SOs calculated using 
model magnitudes and fc-corrections from our BC03:2 model (age 
12 Gyr and r = 1 Gyr). 
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Figure 17. Colour-magnitude relation for E/SOs calculated us- 
ing model magnitudes and fc-corrections the CWW non-evolving 
elliptical template spectrum. 



aa = -(0.1725 ±0.0076). 

Although this evolution rate seems consistent with the A((; — 
r) = —0.17z directly predicted by the model, there is a 
discrepancy in that all the CMR evolution is at 2: < 0.15. 
Indeed, Figures [7] and |5] showed that no single BC03 model 
fits the colours and k-corrections of both the low and high 
redshift E/SOs. 

The result of using the CWW elliptical galaxy template 
to compute the fc-correction is shown in Figure 1171 This 
CMR has 

ao = 0.3312 ± 0.0094 

ai = -(0.01955 ±0.00046) 

02 = -(0.3530 ±0.0076); 

this produces rather more evolution than expected from the 
BC03 models. 

It is evident the CMR and its evolution, as derived 
from the photometry, are extremely sensitive to the adopted 
fc- correct ion. For example, both the BC03 and CWW fc- 
corrections give bluer rest-frame colours than the either the 
spectra or Blanton-Roweis models. 
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Figure 18. CcrR using rest-frame g — r derived from the spectra, 
before and after the signal/noise correction (top and bottom). 
Points corresponding to fewer than 8 galaxies are not plotted. 



5 THE COLOUR-(7 RELATION (CaR) 

The CMR is due to more fundamental correlations between 
colour and velocity dispersion ct, and cr and luminosity 
(Bernardi et al. 2005). The previous section showed that 
colour gradients affect the zero-point and slope of the CMR. 
We now study how the colour-cr relation (hereafter CctR) de- 
pends on the aperture within which the colour is defined. 

Figure [181 shows the CctR for our sample for rest-frame 
g — r from the spectra, before (top) and after (bottom) cor- 
recting for signal/noise. Fitting the galaxies by equation (8), 
but with Mr replaced by log a yields 

ao = 0.3934 ± 0.0041 or 0.3804 ± 0.0040 

ai = 0.18685 ±0.00195 or 0.19061 ± 0.00189 

aa = -(0.1860 ± 0.0049) or - (0.2032 ± 0.0048) 

with rms residuals of 0.0521 magnitudes. 

Note that colour evolution is strongly detected, but it 
is ~ 0.06-0.07z less than in the corresponding CMR. Us- 
ing model g — r and the (uncorrected or corrected) spectra- 
derived k-corrections yields 

ao = 0.3317 ± 0.0043 or 0.3285 ± 0.0044 

ai = 0.19416 ±0.00206 or 0.19139 ± 0.00209 

aa = -(0.2607 ± 0.0052) or - (0.1289 ± 0.0053) 

with rms residuals of 0.0551 magnitudes (see Figure I19|) . 
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Figure 19. Same as previous figure, but now using model g — r 
colour and fc-corrections from the spectra, before and after the 
signal/noise correction. 
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Figure 20. Same as previous figure, but now using fiber colours 
and fc-corrections from the spectra, before and after the sig- 
nal/noise correction. 



This CcrR relation is bluer than that based on the spectra, 
but the slope is almost unchanged. However, it has O.lSz 
less evolution. Finally, using fiber g — r (Figure [20)) instead 
gives 

ao = 0.3740 ± 0.0045 or 0.3709 ± 0.0046 

ai = 0.18186 ±0.00212 or 0.17902 ± 0.00215 

aa = -(0.2859 ± 0.0053) or - (0.1540 ± 0.0040) 

Following Bernardi et al. (2005), if we plot the residual of 
g — r colours (model magnitude and signal/noise corrected) 
to the CMR, {g — r)rf — {{g~r)rf\Mr, z), against the residual 
of log a to the best-fit a — Mr relation, we see the two are 
correlated (Figure I21[) . The correlation coefficient is 0.2994 
with a slope 0.2033 ± 0.0024, an intercept 0.0000 ± 0.0002, 
and a standard deviation 0.0559 mag. However, if we plot the 
residual of colours to the CcrR, [g — r)^/ — ((<? — r)rf\o, z) 
against the residual of M,. to the best-fit Mr — a relation 
(Figure [22])), we find a much lower correlation coefficient 
0.0231 and slope 0.0022 ± 0.0004 (with an intercept again 
0.0000 ±0.0002 and standard deviation 0.0560 magnitudes). 

Again we find that g — r colour in E/SOs are more funda- 
mentally dependent on velocity dispersion than on Mr, and 
furthermore, the absence of an Mr correlation in the second 
case implies that the CcrR is the same for all bins in magni- 
tude at a fixed velocity dispersion. Hence our determination 



of the CcrR in this Section should not be significantly biased 
by our use of a flux-limited sample. 

In both the CMR and CcrR, rest-frame colours from 
the spectra show stronger g — r evolution (0.20-0.27^) than 
model-magnitude colours with the 'corrected' k-correction 
(0.07-0.132:). This is an aperture affect, due to the colour 
gradients which also shift the zero-point of both the CMR 
and CcrR bluewards when model colours are used. However, 
it is notable that, while the CMR based on model colours is 
30% shallower than the CMR based on colours from within 3 
arcsec, this does not happen for the CctR. This suggests that 
the strength of the colour gradient is positively correlated 
with luminosity, but not with velocity dispersion. 

We study this in Figure [23] which shows the difference 
between the rest-frame g — r from fiber and model mag- 
nitudes against the residual of Mr to the best-fit Mr — cr 
relation. We find there is a weak but significant correlation, 
with correlation coefficient 0.069, slope -0.0039 ± 0.0002, 
intercept 0.0111 ± 0.0001, and a standard deviation 0.0327 
magnitudes. The slope d{g—r)/dMr — —0.0039, can account 
for the greater part of the additional steepness of the fiber 
and spectra magnitude CMRs relative to the model magni- 
tude CMR. Similarly, the difference between the rest-frame 
g — r from fiber and model magnitudes can be plotted 
against the residual of log a to the best-fit a — Mr rela- 
tion (Figure I24|l . We note that this residual is positively 
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Figure 21. The residual (g — r)rf — {{g — r)rf\Mr,z) plotted 
against the residual logu — {\oga\Mr,z), as a contour plot with 
the best-fit linear regression. 



Figure 23. The difference in rest-frame g—r colours derived using 
fiber magnitudes and model magnitudes (with k-corrections from 
the spectra), plotted against the residual Mr — {Mr\a,z). 
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Figure 22. The residual (g — r)rf — {{g—r)j.f\'T, z) plotted against 
the residual Mr - {Mr\cr, z). 

correlated to galaxy age (Bernardi et al. 2005; Gallazzi et 
al. 2006). Again we find a correlation, with coefficient 0.098, 
slope -0.0372±0.0014, intercept O.OllldbO.OOOl, and a stan- 
dard deviation 0.0326 magnitudes. These results indicate a 
dependence of colour gradient in E/SOs on the gala^xy prop- 
erties, such as age. This is the subject of a work in progress. 



6 SUMMARY AND DISCUSSION 

1. From the SDSS data release 6, with 367471 galaxy 
spectra, we selected 70378 E/SO (early-type) galaxies at 
Q < z < 0.36, using a combination of both morphological 
(pure de Vaucouleurs profile) and spectral (negative eclass 
from principal component analysis) criteria. We measured 
observer-frame and rest-frame g and r-band magnitudes for 



Figure 24. The diflference in rcst-framc g — r colours derived 
using fiber and model magnitudes (with k-corrections from the 
spectra), plotted against the residual log a — (log a\Mr,z). 



these galaxies by integrating over the SDSS spectra (cor- 
rected for atmospheric and Galactic absorption) weighting 
by the appropriate filter functions. The spectra magnitudes 
of these galaxies were then compared with photometric (3 
arcsec aperture) fiber magnitudes from the SDSS imaging 
data. There weis a large scatter (~ 0.1 mag) between the 
two but on average, the spectra gave systematically redder 
g — r colours than the fiber magnitudes, the difference being 
0.03 0.04 magnitudes at lower redshifts. 

Some, but not all, of this discrepancy is due to seeing: 
the flux in the SDSS spectra is reduced relative to the imag- 
ing photometry by about 0.021 mag per arcsec of seeing in 
g and 0.012 mag in r. As the median seeing is 1.8 arcsec 
this may cause A(g — r) ~ 0.01 mag reddening of the spec- 
tra. There may be a additional small systematic differences 
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in the colour calibrations, consistent with the quoted uncer- 
tainties of up to 0.02 mag (Adelman-Mccarthy et al. 2008). 

2. In the g band, the mean difference between gspec and Qfu, 
decreases significantly (~ 0.05 mag) with increasing redshift, 
although the scatter increases. A similar trend is not seen in 
the r band, potentially compromising estimates of how the 
colours of these objects evolve. The redshift dependence of 
Qspec ~ Qfib with redshift is much greater for galaxies with 
the most negative values of the eclass spectral classifica- 
tion parameter. At z ~ 0.3, such objects have low fiux and 
high noise at the blue end of the g-band, and gspec — gjib 
increases at low signal/noise. We calibrated this dependence 
on eclass and signal/noise, and used it to 'correct' spectra- 
derived magnitudes for these effects. Understanding the un- 
derlying cause for these trends requires further work. For 
example, a comparison of SDSS spectra with best-fitting 
BC03-type models, stacking large numbers together to in- 
crease signal/noise, might help to reveal exactly where in 
the spectra the problem is occurring. 

3. The differences between observer-frame and rest-frame 
magnitudes as measured from the spectra provide estimates 
of fc-corrections in g and r. The k-corrections from the SDSS 
spectra of E/SOs were similar to those from BC03 models 
with ages ~ 12 Gyr, and star-formation exponentially de- 
clining with r = 1 Gyr for the most luminous galaxies and 
r = 2 Gyr for the least luminous. These SF histories and 
their moderate luminosity dependence agree with Jimenez 
et al. (2007). 

4. Rest-frame colours {g — r)rf were measured directly from 
the spectra, and also calculated by subtracting fc-corrections 
from the SDSS fiber or model magnitudes. The first two 
colour estimates are based on fixed angular apertures (about 
3 arcsec); the scale associated with the latter is related to the 
half-light radius, so it varies from object to object. Absolute 
magnitudes Mr were calculated from the fc-corrected model 
magnitudes, and were then corrected for evolution. The 
rest-frame colour-magnitude relation (CMR) and colour-a 
(CctR) relations were measured in a series of redshift inter- 
vals. Both are well-approximated by power-laws that tend 
to shift gradually bluewards with increasing redshift. The 
zero-point, slope and evolution of the CMR associated with 
fixed angular apertures (spectral and fiber colours) differ 
from that found for model magnitudes. We found the slope 
30-40% steeper, and note Mei et al. (2006) and Scodeg- 
gio (2001) found similar differences in the CMRs from fixed 
aperture or from half-light colours. These differences are a 
consequence of radial colour gradients. Gradients complicate 
interpretation of evolution of the CMR from fixed apertures 
- evolution of the CMR based on model magnitudes is easier 
to interpret. 

5. Whereas colour gradients affect the zero-point of the CcrR, 
the slope is much less affected. This suggests that gradients 
are not positively correlated with a - but we find indications 
that they are correlated with residuals in the high luminosity 
and low a direction from the {a\A'Ir,z) relation. 

6. The fe-corrections obtained from fitting BC03 template 
spectra to the observed ugriz fluxes (following Blanton & 
Roweis 2007) depend only slightly on whether one uses 
fiber or model magnitudes (< 0.01 mag at 2 < 0.15). 



7. We also calculated the CMR using k-corrrections derived 
from model or template spectra, rather than the SDSS spec- 
tra. The estimate of rest-frame colour evolution is very sen- 
sitive to the assumed fc-correction, and these models gave 
very different results. The BC03 Model 2 gave a reasonable 
(g — r)r/evolution rate — 0.1662z, but at an uneven rate with 
all the evolution at 2 < 0.15, apparently because a model 
with a single SF history could not simultaneously fit the 
more and less luminous galaxies. The CWW fc-correction, 
which is based on a non-evolving spectrum, gave a CMR 
with strong evolution of —0.35892:, so this did not give self- 
consistent results. The fc-corrections from Blanton & Roweis 
(2007) are smaller than those based on the SDSS spectra, 
and so the CMRs show less evolution - the CMR based 
on model magnitudes showed no evolution. This may be an 
artifact of the fact that the BC03 models assume solar abun- 
dance ratios, and this is unrealistic. 
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Table 1. g-band k-corrections for E/SO galaxies (averaged in redshift intervals), derived from the SDSS spectra (with and without the 
signal/noise correction, the Blanton models (based on fibre or model magnitudes) and the CWW non-evolving elliptical galaxy template. 



Redshift 


Spectra-derived 


Blanton 


model 


CWW elliptical 


range 


uncorrected 


corrected 


fibre mag. 


model mag. 


template 


n n n no 
U.U— U.Uz 


0.0467 


0.0465 


U.U4UO 


n nyi oi 


U.U044 


U.Uz— U.U4 


0.0877 


0.0871 


n Q01 fi 


n n'7'7C; 
U.U/ 10 


n nnno 

u.uyyz 


U.U4— U.UD 


0.1.398 


0.1385 


U. i 1.J0 


n 1 ov/i 

U. (4 


U.loD i 


U.UO— U.Uo 


u . ^uuu 


0.1983 


n 1 n/i /I 


n 1 Q ^n 
U. iooy 


U.zzU 1 


0.08-0.10 


0.2690 


0.2645 


0.2610 


0.2520 


0.2953 


0.10-0.12 


0.3526 


0.3474 


0.3410 


0.3308 


0.3851 


0.12-0.14 


0.4389 


0.4317 


0.4128 


0.4129 


0.4780 


0.14-0.16 


0.5329 


0.5230 


0.5087 


0.5001 


0.5756 


0.16-0.18 


0.6324 


0.6193 


0.6021 


0.5944 


0.6784 


0.18-0.20 


0.7406 


0.7236 


0.6993 


0.6932 


0.7853 


0.20-0.22 


0.8500 


0.8274 


0.7965 


0.7908 


0.8887 


0.22-0.24 


0.9470 


0.9191 


0.8870 


0.8828 


0.9854 


0.24-0.26 


1.0454 


1.0115 


0.9760 


0.9724 


1.0787 


0.26-0.28 


1.1480 


1.1071 


1.0638 


1.0584 


1.1688 


0.28-0.30 


1.2641 


1.2104 


1.1583 


1.1537 


1.2653 


0.30-0.32 


1.3614 


1.3009 


1.2564 


1.2461 


1.3682 


0.32-0.34 


1.4141 


1.3456 


1.3721 


1.3737 


1.4717 


0.34-0.36 


1.6124 


1.5301 


1.4303 


1.4229 


1.5524 


Table 2. As above, for the 


r-band k-corrections 






Redshift 


Spectra-derived 


Blanton 


model 


CWW elliptical 


range 


uncorrected 


corrected 


fibre mag. 


model mag. 


template 


0.0-0.02 


0.0174 


0.0173 


0.0170 


U.Ui / U 


n nonfi 
U.UzUO 
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0.0348 


0.0346 


0.0372 


U.UoDi 


U.Uoo4 


0.04-0.06 


0.0585 


0.0580 


0.0399 


U.Uo/o 


n ^^^^^ 
U.UDUU 


0.06-0.08 


0.0834 


0.0825 


0.0845 


n none; 


n riQOQ 
U.Uozo 


0.08-0.10 


0.1090 


0.1077 


0.1078 


0.1041 


0.1105 


0.10-0.12 


0.1408 


0.1392 


0.1362 


0.1330 


0.1454 


0.12-0.14 


0.1711 


0.1693 


0.1555 


0.1619 


0.1778 


0.14-0.16 


0.1985 


0.1965 


0.5087 


0.5001 


0.5756 


0.16-0.18 


0.2266 


0.2246 


0.2181 


0.2163 


0.2292 


0.18-0.20 


0.2552 


0.2531 


0.2462 


0.2452 


0.2557 


0.20-0.22 


0.2861 


0.2840 


0.2757 


0.2743 


0.2841 


0.22-0.24 


0.3168 


0.3147 


0.3047 


0.3029 


0.3134 


0.24-0.26 


0.3510 


0.3490 


0.3345 


0.3330 


0.3453 


0.26-0.28 


0.3863 


0.3845 


0.3688 


0.3673 


0.3807 


0.28-0.30 


0.4377 


0.4360 


0.4079 


0.4071 


0.4192 


0.30-0.32 


0.4675 


0.4665 


0.4370 


0.4316 


0.4562 


0.32-0.34 


0.5274 


0.5257 


0.4861 


0.4784 


0.4955 


0.34-0.36 


0.5430 


0.5425 


0.5111 


0.4961 


0.5345 



